Initial stages of CO2 adsorption on CaO: a combined experimental and computational study.

Room temperature adsorption of carbon dioxide (CO2) on monocrystalline CaO(001) thin films grown on a Mo(001) substrate was studied by infrared reflection-absorption spectroscopy (IRAS) and quantum chemical calculations. For comparison, CO2 adsorption was examined on poorly ordered, nanoparticulate CaO films prepared on Ru(0001). For both systems, CO2 readily adsorbs on the clean CaO surface. However, additional bands were observable on the CaO/Ru(0001) films compared with CaO/Mo(001), because the stricter IRAS surface selection rules do not apply to adsorption on the disordered thin films grown on Ru(0001). Spectral evolution with increasing exposure of the IRA bands suggested the presence of several adsorption sites which are consecutively populated by CO2. Density functional calculations showed that CO2 adsorption occurs as monodentate surface carbonate (CO32-) species at monatomic step sites and other low-coordinated sites, followed by formation of carbonates on terraces, which dominate at increasing CO2 exposure. To explain the coverage-dependent IRAS results, we propose CO2 surface islanding from the onset, most likely in the form of pairs and other chain-like species, which were calculated as thermodynamically favorable. The calculated adsorption energy for isolated CO2 on the terrace sites (184 ± 10 kJ mol-1) is larger than the adsorption energy obtained by temperature programmed desorption (∼120-140 kJ mol-1) and heat of adsorption taken from microcalorimetry measurements at low coverage (∼125 kJ mol-1). However, the calculated adsorption energies become less favorable when carbonate chains intersect on CaO terraces, forming kinks. Furthermore, our assignments of the initial stages of CO2 adsorption are consistent with the observed coverage effect on the CO2 adsorption energy measured by microcalorimetry and the IRAS results.

Computational Study of Fluorinated Diglyoxime-Iron Complexes: Tuning the Electrocatalytic Pathways for Hydrogen Evolution.

The ability to tune the properties of hydrogen-evolving molecular electrocatalysts is important for developing alternative energy sources. Fluorinated diglyoxime-iron complexes have been shown to evolve hydrogen at moderate overpotentials. Herein two such complexes, [(dAr(F)gBF2)2Fe(py)2], denoted A, and [(dAr(F)g2H-BF2)Fe(py)2], denoted B [dAr(F)g = bis(pentafluorophenyl-glyoximato); py = pyridine], are investigated with density functional theory calculations. B differs from A in that one BF2 bridge is replaced by a proton bridge of the form O-H-O. According to the calculations, the catalytic pathway for A involves two consecutive reduction steps, followed by protonation of an Fe(0) species to generate the active Fe(II)-hydride species. B is found to proceed via two parallel pathways, where one pathway is similar to that for A, and the additional pathway arises from protonation of the O-H-O bridge, followed by spontaneous reduction to an Fe(0) intermediate and intramolecular proton transfer from the ligand to the metal center or protonation by external acid to form the same active Fe(II)-hydride species. Simulated cyclic voltammograms (CVs) based on these mechanisms are in qualitative agreement with experimental CVs. The two parallel pathways identified for B arise from an equilibrium between the protonated and unprotonated ligand and result in two catalytic peaks in the CVs. The calculations predict that the relative probabilities for the two pathways, and therefore the relative magnitudes of the catalytic peaks, could be tuned by altering the pK(a) of the acid or the substituents on the ligands of the electrocatalyst. The ability to control the catalytic pathways through acid strength or ligand substituents is critical for designing more effective catalysts for energy conversion processes.

Nickel phlorin intermediate formed by proton-coupled electron transfer in hydrogen evolution mechanism.

The development of more effective energy conversion processes is critical for global energy sustainability. The design of molecular electrocatalysts for the hydrogen evolution reaction is an important component of these efforts. Proton-coupled electron transfer (PCET) reactions, in which electron transfer is coupled to proton transfer, play an important role in these processes and can be enhanced by incorporating proton relays into the molecular electrocatalysts. Herein nickel porphyrin electrocatalysts with and without an internal proton relay are investigated to elucidate the hydrogen evolution mechanisms and thereby enable the design of more effective catalysts. Density functional theory calculations indicate that electrochemical reduction leads to dearomatization of the porphyrin conjugated system, thereby favoring protonation at the meso carbon of the porphyrin ring to produce a phlorin intermediate. A key step in the proposed mechanisms is a thermodynamically favorable PCET reaction composed of intramolecular electron transfer from the nickel to the porphyrin and proton transfer from a carboxylic acid hanging group or an external acid to the meso carbon of the porphyrin. The C-H bond of the active phlorin acts similarly to the more traditional metal-hydride by reacting with acid to produce H2. Support for the theoretically predicted mechanism is provided by the agreement between simulated and experimental cyclic voltammograms in weak and strong acid and by the detection of a phlorin intermediate through spectroelectrochemical measurements. These results suggest that phlorin species have the potential to perform unique chemistry that could prove useful in designing more effective electrocatalysts.

Role of pendant proton relays and proton-coupled electron transfer on the hydrogen evolution reaction by nickel hangman porphyrins.

The hangman motif provides mechanistic insights into the role of pendant proton relays in governing proton-coupled electron transfer (PCET) involved in the hydrogen evolution reaction (HER). We now show improved HER activity of Ni compared with Co hangman porphyrins. Cyclic voltammogram data and simulations, together with computational studies using density functional theory, implicate a shift in electrokinetic zone between Co and Ni hangman porphyrins due to a change in the PCET mechanism. Unlike the Co hangman porphyrin, the Ni hangman porphyrin does not require reduction to the formally metal(0) species before protonation by weak acids in acetonitrile. We conclude that protonation likely occurs at the Ni(I) state followed by reduction, in a stepwise proton transfer-electron transfer pathway. Spectroelectrochemical and computational studies reveal that upon reduction of the Ni(II) compound, the first electron is transferred to a metal-based orbital, whereas the second electron is transferred to a molecular orbital on the porphyrin ring.

Proton-coupled electron transfer in molecular electrocatalysis: theoretical methods and design principles.

Molecular electrocatalysts play an essential role in a wide range of energy conversion processes. The objective of electrocatalyst design is to maximize the turnover frequency and minimize the overpotential for the overall catalytic cycle. Typically, the catalytic cycle is dominated by key proton-coupled electron transfer (PCET) processes comprised of sequential or concerted electron and proton transfer steps. Theoretical methods have been developed to investigate the mechanisms, thermodynamics, and kinetics of PCET processes in electrocatalytic cycles. Electronic structure methods can be used to calculate the reduction potentials and pKa's and to generate thermodynamic schemes, free energy reaction pathways, and Pourbaix diagrams, which indicate the most stable species under certain conditions. These types of calculations have assisted in identifying the thermodynamically favorable mechanisms under specified experimental conditions, such as acid strength and overpotential. Such calculations have also revealed linear correlations among the thermodynamic properties, which can be used to predict the impact of modifying the ligands, substituents, or metal centers. The thermodynamic properties can be tuned with electron-withdrawing or electron-donating substituents. Ligand modification can exploit the role of noninnocent ligands. For example, ligand protonation typically decreases the overpotential. Calculations of rate constants for electron and proton transfer, as well as concerted PCET, have assisted in identifying the kinetically favorable mechanisms under specified conditions. The concerted PCET mechanism is thought to lower the overpotential required for catalysis by avoiding high-energy intermediates. Rate constant calculations have revealed that the concerted mechanism involving intramolecular proton transfer will be favored by designing more flexible ligands that facilitate the proton donor-acceptor motion while also maintaining a sufficiently short equilibrium proton donor-acceptor distance. Overall, theoretical methods have assisted in the interpretation of experimental data and the design of more effective molecular electrocatalysts.

Effects of ligand modification and protonation on metal oxime hydrogen evolution electrocatalysts.

The design of hydrogen-evolving electrocatalysts that operate at modest overpotentials is important for solar energy devices. The M(II/I) reduction potential for metal diimine-dioxime and diglyoxime electrocatalysts is often related to the overpotential required for hydrogen evolution. Herein the impact of ligand modification and protonation on the M(II/I) reduction potentials for cobalt, nickel, and iron diimine-dioxime and diglyoxime complexes is investigated with computational methods. The calculations are consistent with experimental data available for some of these complexes and additionally provide predictions for complexes that have not yet been synthesized. The calculated pKa's imply that ligand protonation is likely to occur at the O-H-O bridge but not at other ligand sites for these complexes. Moreover, the calculations imply that a ligand-protonated Co(III)-hydride intermediate is formed along the H2 production pathway for catalysts containing an O-H-O bridge in the presence of sufficiently strong acid. The calculated M(II/I) reduction potentials indicate that the anodic shift due to protonation of the O-H-O bridge is greater than that due to replacing the O-H-O bridge with an O-BF2-O bridge for cobalt and nickel but not for iron complexes. Experiments suggest degradation for complexes with two O-H-O bridges and alternative mechanisms for certain iron complexes with two O-BF2-O bridges. Asymmetric cobalt, nickel, and strongly electron withdrawing substituted iron diimine-dioxime and diglyoxime complexes containing a single O-H-O bridge are proposed to be effective hydrogen evolution electrocatalysts with relatively low overpotentials in acetonitrile and water. These insights are important for the design of efficient aqueous-based hydrogen-evolving catalysts.

Photoinduced proton-coupled electron transfer of hydrogen-bonded p-nitrophenylphenol-methylamine complex in solution.

Proton-coupled electron transfer can occur through concerted (electron-proton transfer, EPT) or sequential mechanisms, but this distinction becomes less well-defined for photoinduced reactions. These issues have been examined with transient absorption experiments on a hydrogen-bonded complex consisting of p-nitrophenylphenol and t-butylamine. These experiments revealed two spectroscopically distinct states: the higher-energy excited state was interpreted to be a conventional intramolecular charge transfer (ICT) state within the p-nitrophenylphenol, whereas the lower-energy state was interpreted to be an ICT-EPT state, where photoexcitation resulted in both ICT and the shifting of electronic density corresponding to effective proton transfer from the phenol to the amine. In the present work, the singlet excited states of the hydrogen-bonded p-nitrophenylphenol-methylamine complex in 1,2-dichloroethane are studied with time-dependent density functional theory and higher-level ab initio methods. The calculations suggest that the ππ* state, which is the S(1) state at the Franck-Condon geometry, corresponds to the state denoted ICT-EPT in the experimental analysis, whereas the nπ* state, which is the S(2) state at this geometry, likely corresponds to the state denoted ICT in the experimental analysis. According to the calculations, the ππ* state has charge-transfer character, as well as a change in electronic density on the amine, with the minimum-energy structure corresponding to the proton bonded to the nitrogen acceptor, consistent with proton transfer. The nπ* state has little charge-transfer character, as well as negligible change in electronic density on the amine, with the minimum-energy structure corresponding to the proton bonded to the oxygen donor. The calculations also provide evidence of an avoided crossing between these two states located energetically close to the Franck-Condon point. These calculations provide the foundation for future nonadiabatic molecular dynamics studies of the relaxation process.

Computational study of anomalous reduction potentials for hydrogen evolution catalyzed by cobalt dithiolene complexes.

The design of efficient hydrogen-evolving catalysts based on earth-abundant materials is important for developing alternative renewable energy sources. A series of four hydrogen-evolving cobalt dithiolene complexes in acetonitrile-water solvent is studied with computational methods. Co(mnt)(2) (mnt = maleonitrile-2,3-dithiolate) has been shown experimentally to be the least active electrocatalyst (i.e., to produce H(2) at the most negative potential) in this series, even though it has the most strongly electron-withdrawing substituents and the least negative Co(III/II) reduction potential. The calculations provide an explanation for this anomalous behavior in terms of protonation of the sulfur atoms on the dithiolene ligands after the initial Co(III/II) reduction. One fewer sulfur atom is protonated in the Co(II)(mnt)(2) complex than in the other three complexes in the series. As a result, the subsequent Co(II/I) reduction step occurs at the most negative potential for Co(mnt)(2). According to the proposed mechanism, the resulting Co(I) complex undergoes intramolecular proton transfer to form a catalytically active Co(III)-hydride that can further react to produce H(2). Understanding the impact of ligand protonation on electrocatalytic activity is important for designing more effective electrocatalysts for solar devices.

Substituent effects on cobalt diglyoxime catalysts for hydrogen evolution.

The design of efficient, robust, and inexpensive hydrogen evolution catalysts is important for the development of renewable energy sources such as solar cells. Cobalt diglyoxime complexes, Co(dRgBF(2))(2) with substituents R, are promising candidates for such electrocatalysts. The mechanism for hydrogen production requires a series of reduction and protonation steps for various monometallic and bimetallic pathways. In this work, the reduction potentials and pK(a) values associated with the individual steps were calculated for a series of substituents. The calculations revealed a linear relation between the reduction potentials and pK(a) values with respect to the Hammett constants, which quantify the electron donating or withdrawing character of the substituents. Additionally, the reduction potentials and pK(a) values are linearly correlated with each other. These linear correlations enable the prediction of reduction potentials and pK(a) values, and thus the free energy changes along the reaction pathways, to assist in the design of more effective cobaloxime catalysts.

Theoretical analysis of mechanistic pathways for hydrogen evolution catalyzed by cobaloximes.

The mechanistic pathways for hydrogen evolution catalyzed by cobalt complexes with supporting diglyoxime ligands are analyzed with computational methods. The cobaloximes studied are Co(dmgBF(2))(2) (dmg = dimethylglyoxime) and Co(dpgBF(2))(2) (dpg = diphenylglyoxime) in acetonitrile. The reduction potentials and pK(a) values are calculated with density functional theory in conjunction with isodesmic reactions, incorporating the possibility of axial solvent ligand loss during the reduction process. The solvent reorganization energies for electron transfer between the cobalt complex and a metal electrode and the inner-sphere reorganization energies accounting for intramolecular rearrangements and the possibility of ligand loss are also calculated. The relative reduction potentials agree quantitatively with the available experimental values. The pK(a)s and reorganization energies agree qualitatively with estimates based on experimental data. The calculations suggest that a peak measured at ca. -1.0 V vs SCE in cyclic voltammetry experiments for Co(dmgBF(2))(2) is more likely to correspond to the Co(II)H/Co(I)H reduction potential than the Co(III)H/Co(II)H reduction potential. The calculations also predict pK(a) values of Co-hydride complexes and reduction potentials for both cobaloximes that have not been determined experimentally. The results are consistent with a mechanism in which the Co(III) and Co(II) complexes have two axial solvent ligands and the Co(I) complex has a single axial ligand along the reaction pathway. Analysis of the free energy diagrams generated for six different monometallic and bimetallic hydrogen production pathways identified the most favorable pathways for Co(dmgBF(2))(2) and tosic acid. The thermodynamically favored monometallic pathway passes through a Co(III)H intermediate, and Co(II)H reacts with the acid to produce H(2). The thermodynamically favored bimetallic pathways also pass through the Co(III)H intermediate, but the pathways in which two Co(III)H or two Co(II)H complexes react to produce H(2) are not thermodynamically distinguishable with these methods. On the basis of the electrostatic work term associated with bringing the two cobalt complexes together in solution, the preferred bimetallic pathway involves the reaction of two Co(III)H complexes to produce H(2). This mechanistic insight is important for designing more effective catalysts for solar energy conversion.